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Abstract. The exclusive electroproduction of tt + above the resonance region was studied using the CEBAF 
Large Acceptance Spectrometer (CLAS) at Jefferson Laboratory by scattering a 6 GeV continuous electron 
beam off a hydrogen target. The large acceptance and good resolution of CLAS, together with the high 
luminosity, allowed us to measure the cross section for the 7*p — > niv + process in 140 (Q 2 , xb, i) bins: 
0.16 < x B < 0.58, 1.6 GeV 2 < Q 2 < 4.5 GeV 2 and 0.1 GeV 2 < -t < 5.3 GeV 2 . For most bins, the statistical 
accuracy is on the order of a few percent. Differential cross sections are compared to four theoretical 
models, based either on hadronic or on partonic degrees of freedom. The four models can describe the 
gross features of the data reasonably well, but differ strongly in their ingredients. In particular, the model 
based on Generalized Parton Distributions (GPDs) contain the interesting potential to experimentally 
access transversity GPDs. 
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1 Introduction 

One of the major challenges in contemporary nuclear physics 
is the study of the transition between hadronic and par- 
tonic pictures of the strong interaction. At asymptoti- 
cally short distances, the strong force is actually weak 
and the appropriate degrees of freedom are the quarks 
and gluons (partons) whose interaction can be quantified 
very precisely by perturbative Quantum Chromodynamics 
(pQCD). However, at larger distances on the order of one 
Fermi, effective theories that take hadrons as elementary 
particles whose interactions are described by the exchange 
of mesons appear more applicable. The connection be- 
tween these two domains is not well understood. In order 
to make progress, a systematic study of a series of hadronic 
reactions probing these intermediate distance scales is nec- 
essary. The exclusive electroproduction of a meson (or of 
a photon) from a nucleon, j*N —> N'M, is particularly 
interesting. Indeed, it offers two ways to vary the scale 
of the interaction and therefore to study this transition 
regime. One can vary the virtuality of the incoming pho- 



ton Q 2 — — (p e — p' e ) 2 , which effectively represents the 
transverse size of the probe, or the momentum transfer to 
the nucleon t — (pn — p'n) 2 > which effectively represents 
the transverse size of the target. Here, p e and p' e are the 
initial and scattered electron four- momenta and pn and 
p' N are the initial and final nucleon four-momenta, respec- 
tively. Figure Q] sketches the transition regions that have 
been experimentally explored until now (lightly shaded 
areas) as a function of these two variables, Q 2 and \t\. In 
this figure, we keep, quite arbitrarily, only the experiments 
for which |i| > 3 GeV 2 in photoproduction (SLAC [Tj and 
JLab [2]) and Q 2 > 1.5 GeV 2 in electroproduction (Cor- 
nell (H], JLab [SPlfT] and HERMES 0). These are the 
domains for which, we believe, there are chances to observe 
first signs that partonic degrees of freedom play a role in 
the reactions. The darkly shaded area in Fig. [1] represents 
the phase space covered by the present work. It is divided 
into 140 (Q 2 , xb or W, t) bins, to be compared to only a 
few (Q 2 , xb or W, t) bins in the lightly shaded areas for 
the previous electroproduction experiments. 
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We also display in Fig.[T]three Feynman-type diagrams 
illustrating the mechanisms believed to be at stake for the 
7*7V — > N'M process: at asymptotically large-Q 2 , asymp- 
totically large- \t | (both in terms of partonic degrees of 
freedom) and at low-Q 2 and low-|t| (in terms of hadronic 
degrees of freedom). 

At asymptotically large Q 2 and small \t\ (along the 
vertical axis in Fig. [1]), the exclusive electroproduction of 
a meson should be dominated by the so-called "handbag 
diagram" [9lllOUllU12j . The initial virtual photon hits a 
quark in the nucleon and this same quark, after a sin- 
gle gluon exchange, ends up in the final meson. A QCD 
factorization theorem [12] states that the complex quark 
and gluon non-perturbative structure of the nucleon is de- 
scribed by the Generalized Parton Distributions (GPDs). 
For the ir + channel at leading twist in QCD, i.e. at asymp- 
totically large Q 2 , the longitudinal part of the cross sec- 
tion a i, is predicted to be dominant over the transverse 
part <jt- Precisely, duhldt should scale as l/Q 6 at fixed 
xb and \t\, while dar/dt should scale as l/Q 8 . It is pre- 
dicted that <7l is sensitive to the helicity-dependent GPDs 
E and H [T2] while, if higher-twist effects are taken into 
account and factorization is assumed, ut is sensitive to 
the transversity GPDs, H T and E T = 2H T + Et [15]. 



CM 
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Fig. 1. (color online). Schematic representation of the 7*iV — > 
N'n process (above the resonance region) in different regions 
of the (Q 2 , t) plane. The Feynman diagrams describe the re- 
action in terms of meson exchanges at low Q 2 and \t\, in terms 
of GPDs at large Q 2 and small \t\, and in terms of hadron 
distribution amplitudes (DA) at large \t\. The lightly shaded 
areas (magenta and green online) represent approximately the 
experimentally explored regions up to now. The darkly shaded 
area (blue online) represents the phase space covered by this 
work. 



At large values of \t\, in photoproduction (i.e. along the 
horizontal axis in Fig. [T]) but also presumably in electro- 
production, the 7^*^iV — > N'M process should be domi- 
nated by the coupling of the (virtual) photon to one of the 
valence quarks of the nucleon (or of the produced meson) , 
with minimal interactions among the valence quarks. In 
this regime, a QCD factorization theorem states that the 
complex structure of the hadrons is described by distribu- 
tion amplitudes (DA) which at small distances (large \t\) 
can be reduced to the lowest Fock states, i.e. 3 quarks for 
the nucleon and q-q for the meson |14) . At sufficiently high 
energy, constituent counting rules (CCR) [15] predict an 
s~ 7 scaling of the differential cross section da/dt at fixed 
center-of-mass pion angles, provided |s|, \t\, and \u\ are 
all large. Here s = W 2 is the squared invariant mass of 
the 7*-p system and u = (p* — p' N ) 2 is given in terms of 
the four-vectors p* — p e — p' e and p' N for the final-state 
nucleon. The large \t\ and |w| region corresponds typically 
to a center-of-mass pion angle 6 cm ps 90°. In this domain, 
the CCR predict da/dt = f(0 cm )s 2 ~ n for the energy de- 
pendence of the cross section, where f(0 cm ) depends on 
details of the dynamics of the process and n is the total 
number of point-like particles and gauge fields in the ini- 
tial and final states. For example, our reaction j*p — > nir + 
should have n — 9, since there is one initial photon, three 
quarks in the initial and the final nucleons, and two in the 
final pion. 

Many questions are open, in particular at which Q 2 
and s do such scaling laws start to appear. Even if these 
respective scaling regimes are not reached at the present 
experimentally accessible Q 2 and s values, can one never- 
theless extract GPDs or DAs, provided that some correc- 
tions to the QCD leading-twist mechanisms are applied? 
Only experimental data can help answer such questions. 



2 Insights from previous experiments with 
respect to partonic approaches 

The two most recent series of experiments that have mea- 
sured exclusive ir + electroproduction off the proton, in 
the large-Q 2 , low-|f| regime where the GPD formalism is 
potentially applicable, have been conducted in Hall C at 
Jefferson Lab (JLab) [SH7] and at HERMES [5]. 

The Hall C experiments, with 2 to 6 GeV electron 
beam energies, separated the ox and <jt cross sections 
of the 7*p — > 7i7T + process using the Rosenbluth technique 
for 0.17 < x B < 0.48 and Q 2 up to 3.91 GeV 2 . The term 
<tl dominated the cross section for \t\ < 0.2 GeV 2 , while 
cjt was dominant for larger \t\ values. These data were 
compared to two GPD-based calculations, hereafter re- 
ferred to as VGG 16] and GK (BET] from the initials of 
the models' authors. The comparison of the data with the 
VGG model can be found in the Hall C publications [5,6 
while the comparison with the GK model can be found in 
the GK publications [T3"Ill7| . For <tl, which should be the 
QCD leading-twist contribution, these GPD calculations 
were found to be in general agreement with the magni- 
tude and the Q 2 - and t- dependencies of the experimental 
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data. In these two calculations the main contribution to 
stems from the E GPD, which is modeled either entirely 
as pion-exchange in the i-channel [16j or is at least dom- 
inated by it [T5UT?] (see Refs. [T51IT5] for the connection 
between the i-channel pion-exchange and the E GPD). 
This term is also called the "pion pole", and the differ- 
ence between the two calculations lies in the particular 
choice made for the i-channel pion propagator (Reggeized 
or not) and the introduction of a hadronic form factor 
or not at the wNN vertex. In both calculations, ctl con- 
tains higher-twist effects because the pure leading-twist 
component of the pion pole largely underestimates the 
data. Only the GK model, which explicitly takes into ac- 
count higher-twist quark transverse momentum, is able 
to calculate ot- Agreement between data and calculation 
is found only if the Ht transversity GPD is introduced, 
which makes up most of <tt- 

The HERMES experiment used 27.6 GeV electron and 
positron beams to measure the 7*p — > mt + cross section 
at four (xb, Q 2 ) values, with xb ranging from 0.08 to 0.35 
and Q 2 from 1.5 to 5 GeV 2 . Since all data were taken at a 
single beam energy, no longitudinal/transverse separation 
could be carried out. The differential cross section dajdt 
was compared to the same two GPD models mentioned 
above. The GK model, which calculates both the longitu- 
dinal and transverse parts of the cross section, displays the 
same feature as for the lower energy JLab data, i.e. a dom- 
inance of Oh up to — i « 0.2 GeV 2 , after which ar takes 
over. The sum of the transverse and longitudinal parts of 
the cross section calculated by the GK model is in very 
good agreement with the data over most of the i range 
measured at HERMES \TM2\- The VGG model, which 
calculates only the longitudinal part of the cross section, 
is in agreement with the data only for low t values [5]. 
Again, in both calculations, ol is dominated by the E 
GPD, modeled essentially by the pion pole term, and or, 
in the GK model, is due to the transversity GPDs. The 
HERMES experiment also measured the transverse target 
spin asymmetry Ajjt for the 7*p — » mr + process, which 
indicate [131117] that the transversity GPDs Ht or Et in- 
deed play an important role in the process, confirming the 
approach of the GK group. 

The comparison between the JLab Hall C and HER- 
MES experiments and the two GPD-based calculations 
yields very encouraging signs that, although higher- twist 
contributions definitely play a major role, these data can 
be interpreted in terms of GPDs, in particular transversity 
GPDs. More precise and extensive data would be highly 
useful to confirm these findings. Firstly, the present CLAS 
experiment extends somewhat the (xb, Q 2 ) phase space 
previously covered by the JLab Hall C experiments and 
secondly, it covers 20 (xb, Q 2 ) bins (with statistical er- 
rors of a few percent on average) which doubles the num- 
ber of bins of the JLab Hall C experiments (and triples 
the HERMES number of bins). These new data are im- 
portant to test the present GPD-based model calculations 
and, if successful, bring more stringent constraints on the 
current GPD parametrizations. 



The large- |i | (large- domain, where the DA formal- 
ism is asymptotically applicable for 7 ( -*-*p — » nit + , has 
so far been explored only in high-energy photoproduction 
at SLAC [T) and intermediate-energy photoproduction at 
JLab [10] ■ While the SLAC data tend to follow the s~ 7 
scaling asymptotic prediction, for a 90° center-of-mass an- 
gle, the more recent JLab data, which are compatible with 
the SLAC data but are more precise, actually reveal some 
large oscillations around this s~ 7 behavior. 

In recent years a similar trend, i.e. "global" scaling be- 
havior, has been observed in deuteron photo-disintegration 
experiments 21,22,231124]. and also in hyperon photopro- 
duction [25] . It would be interesting to see this in exclusive 
pion electroproduction and if so, whether the oscillations 
disappear as Q 2 increases. The measurement presented in 
this article is the first one to explore this large- \t\, large- 
ly domain (6 cm ps 90°) for y/s > 2 GeV in ir + exclusive 
electroproduction off the proton. The present CLAS elec- 
troproduction experiment covers a i-range up to ~ 5 GeV 2 
while the largest |i|-values measured by Hall C are « 0.9 
GeV 2 and by HERMES 2 GeV 2 . 

3 The Experiment 




Fig. 2. Three-dimensional view of the CLAS detector. 



The measurement was carried out with the CEBAF 
Large Acceptance Spectrometer (CLAS) 26 . A schematic 
view of CLAS is shown in Fig. [5] CLAS has a toroidal 
magnetic field generated by six flat superconducting coils 
(main torus), arranged symmetrically around the azimuth. 
Six identical sectors are independently instrumented with 
34 layers of drift cells for particle tracking (Rl, R2, R3), 
plastic scintillation counters for time-of-flight (TOF) mea- 
surements, gas threshold Cherenkov counters (CC) for 
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electron and pion separation, and electromagnetic calorime- 
ters (EC) for photon and neutron detection. To aid in elec- 
tron/pion separation, the EC is segmented into an inner 
part closer to the target and an outer part further away 
from the target. CLAS covers on average 80% of the full 
47T solid angle for the detection of charged particles. The 
azimuthal acceptance is maximum at a polar angle of 90° 
and decreases at forward angles. The polar angle coverage 
ranges from about 8° to 140° for the detection of tt + . The 
scattered electrons are detected in the CC and EC, which 
extend from 8° to 45°. 

The target is surrounded by a small toroidal magnet 
(mini-torus). This magnet is used to shield the drift cham- 
bers closest to the target from the intense low-energy elec- 
tron background resulting from M0ller scattering. 

A Faraday cup, composed of 4000 kg of lead and 75 
radiation lengths thick, is located in the beam dump, 29 
meters downstream the CLAS target. It completely stops 
the electrons and thus allows to measure the accumulated 
charge of the incident beam and therefore the total flux 
of the beam . 

The specific experimental data set "el-6" used for this 
analysis was collected in 2001. The incident beam had an 
average intensity of 7 nA and an energy of 5.754 GeV. 
The 5-cm-long liquid-hydrogen target was located 4 cm 
upstream of the CLAS center. This offset of the target 
position was found to optimize the acceptance of forward- 
going positively charged particles. The main torus magnet 
was set to 90% of its maximum field, corresponding to 
an integral magnetic field of « 2.2 Tm in the forward 
direction. The torus current during the run was very stable 
(< 0.03%). Empty-target runs were performed to measure 
contributions from the target cell windows. 

In this analysis, the scattered electron and the pro- 
duced 7r + were detected and the final state neutron deter- 
mined from missing mass. The continuous electron beam 
provided by CEBAF is well suited for measurements in- 
volving two or more final-state particles in coincidence, 
leading to very small accidental coincidence contributions, 
smaller than 10~ 3 , for the instantaneous luminosity of 
10 cm _2 s _1 of the present measurement. 

Raw data were subjected to the calibration and recon- 
struction procedures that are part of the standard CLAS 
data analysis sequence. Stringent kinematic cuts were ap- 
plied to select events with one electron candidate and only 
one positively charged track. These events were then sub- 
jected to further selection criteria described in the follow- 
ing Section. Throughout the analysis, the experimental 
data distributions were compared to the output of our 
Monte Carlo program GSIM (see Sec.!]). 

A schematic illustration of electron scattering off a nu- 
clcon target producing an outgoing nucleon and one pion 
is shown in Fig. [3] The scattered electron angle 9 e is given 
in the laboratory frame. The angle between the virtual 
photon three-momentum and the direction of the pion is 
denoted as 9% and the angle between the electron scat- 
tering plane and hadronic production plane is denoted as 
0*. These two angles are defined in the center-of-mass 
frame of the hadronic system. The angle 0* is defined so 




SCATTERING PLANE 



Fig. 3. Kinematics of exclusive single tv^ 
from a proton target. 



electroproduction 



that the scattered electron lies in the </>* = 0° half plane 
with the z-axis pointing along the virtual photon momen- 
tum. For exclusive single n + production from the proton, 
we request the simultaneous detection of one single elec- 
tron and of one single ir + in CLAS and the final state 
neutron will be identified by the missing mass squared 
((p e +Pn) — (p' e +Ptt)) 2 , where p^ is the four- momentum 
of the detected 7r + . The kinematic range and bin sizes are 
chosen to provide reasonable statistics in each bin. These 
are summarized in Table [1] 



Table 1. Kinematic bins used in this analysis. 



Variable Number of bins 



Range 



Bin size 



■I'D 



Q 2 



0.16 - 0.58 



1.6 

3.1 



3.1 GeV 
4.5 GeV 2 



0.1 - 1.9 GeV 2 
1.9 - 4.3 GeV 2 
4.3 - 5.3 GeV 2 



0.06 

0.3 GeV 2 
0.5 GeV 2 

0.3 GeV 2 
0.8 GeV 2 
1.0 GeV 2 



Our aim is to extract the three-fold differential cross 



section 



r dQ 2 dx B dt 



where: 



1 



d 3 <7 



n w (Q 2 ,x B ,t) 



rdQ 2 dx B dt C mt AQ 2 Ax B At 1 



(1) 



with 



n w (Q 2 ,XB,t) is the weighted number of ep —> e'mr + 
events in a given bin (Q 2 , x B , t); in particular, n w {x B , 
Q 2 , —t) contains the detector's acceptance correction 
factor Acc(x B , Q 2 , —t, <j>*) (see Sec. I5.1[) and the cor- 
rection factor due to radiative effects F ra d{x B , Q 2 , —t) 
(see Sec.m]), 

Li n t is the effective integrated luminosity, 
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— AQ 2 . Axb and At are the corresponding bin widths 
(see Table [I]); for bins not completely hlled, because of 
W or E' cuts on the electron for instance (see fig. H]), 
the phase space AQ 2 AxsAt includes a volume cor- 
rection and the Q 2 and xb central values are modified 
accordingly. 

In the following three sections, we detail the various cuts 
and correcting factors entering the definition of n w (Q 2 , xb, 



4 Data Analysis 

4.1 Particle identification and event selection 

4.1.1 Electron identification 

The electrons are identified at the trigger level by requir- 
ing at least 640 MeV energy deposited in the EC in co- 
incidence with a signal in the CC (which triggers on one 
photoelectron). 

Additional requirements for particle identification (PID) 
were used in the off-line analysis to refine the electron 
identification. First, we required that the EC and CC hits 
matched with a reconstructed track in the drift cham- 
bers (DC). Second, we correlated the energy deposited in 
the EC and the momentum obtained by the track recon- 
struction in the DC. This is aimed at removing the pion 
contamination. Electrons deposit energy in proportion to 
their incident energy in the calorimeter whereas pions are 
minimum ionizing and deposit a constant amount of en- 
ergy in the calorimeter. The ratio of the total deposited 
energy in the EC to the momentum of the particle is called 
the sampling fraction. For electrons, approximately 30% 
of the total energy deposited in the EC is directly mea- 
sured in the active scintillator material. The remainder 
of the energy is deposited in the lead sheets interleaved 
between the scintillators. Figure 2] shows the sampling 
fraction E/p e versus particle momentum p e . The average 
sampling fraction for electrons was found to be 0.291 for 
this experiment. The solid lines in Fig. @] show the ±3er 
sampling fraction cuts used in this analysis. 

To further reject pions, we required the energy de- 
posited in the inner EC to be larger than 50 MeV. Min- 
imum ionizing particles lose less than this amount in the 
15 cm thickness of the inner EC. 

Fiducial cuts were applied to exclude the EC detector 
edges. When an electron hit is close to an edge, part of the 
shower leaks outside the device; in this case, the energy 
cannot be fully reconstructed from the calorimeter infor- 
mation alone. This problem can be avoided by selecting 
only those electrons lying inside a fiducial volume within 
the EC that excludes the edges. A GEANT-based sim- 
ulation (GSIM) was used to determine the EC-response 
with full electron energy reconstruction. The calorimeter 
fiducial volume was defined by cuts that excluded the in- 
efficient detector regions. 

Particle tracks were reconstructed using the drift cham- 
ber information, and each event was extrapolated to the 
target center to obtain a vertex location. We demanded 




P e [GeV] 

Fig. 4. (color online). EC sampling fraction versus particle 
momentum for the experimental data before (top) and after 
(bottom) EC energy cuts. The solid curves show the ±3cr sam- 
pling fraction cuts which are applied to select electrons. 



that the reconstructed z-vertex position (distance along 
the beam axis from the center of CLAS, with negative val- 
ues indicating upstream of the CLAS center) lies in the 
range —80 mm < z vtx < —8 mm. This is slightly larger 
than the target cell size in order to take into account the 
resolution effects on the vertex reconstruction. 

Finally, a lower limit on the number of photoelectrons 
detected in the photomultiplier tubes of the CC provided 
an additional cut to improve electron identification. The 
number of photoelectrons detected in the CC follows a 
Poisson distribution modified for irregularities in light col- 
lection efficiency for the individual elements of the array. 
For this experiment, a good electron event was required 
to have 3 or more photoelectrons detected in the CC. The 
efficiency of the CC cut was determined from the exper- 
imental data. We fit the number of photoelectrons using 
the modified Poisson distribution. The efficiency range af- 
ter the CC cut is 78% to 99% depending on the kinematic 
region. The correction is then the integral below the cut 
divided by the total integral of the resulting fit function. 
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4.1.2 Positively charged pion identification 

The main cuts to select the 7r + are based on charge, z- 
vertex, fiducial cuts and velocity versus momentum cor- 
relations. The velocity j3 is calculated from the ratio of 
the path length of the reconstructed track, to the time of 
flight. 

Figure [5] shows the /3 versus p distribution for posi- 
tively charged particles from experimental data (top) and 
from the GSIM Monte Carlo simulation (bottom) . A Gaus- 
sian is fit to (3 for bins in momentum p n . A ±1.5(7 cut on 
P is chosen for pion candidates as shown in Fig. [5] (solid 
curves in the plot). Pions and positrons (f3 — 1) are well 
separated below p^ = 250 MeV/c of momentum in the 
experimental data, but this is no longer the case at mo- 
menta larger than 400 MeV/c. For this reason, positrons 
can be mis-identified as pions, which increases the back- 
ground. At higher momenta, there can also be some par- 
ticle mis-identification from protons and kaons. We es- 
timated that the missing mass and vertex cuts reduce 
this mis-identification to the 5 - 10% level. This residual 
background contamination was subtracted as described in 
Secil 

4.2 Fiducial cuts 

4.2.1 Electron fiducial cuts 

The fiducial cuts for electrons were developed to exclude 
regions with non-uniform detector efficiency such as the 
edges of a sector in the CC and EC. The fiducial cut is 
a function of the angles 9 e , <f> e , and momentum p e of the 
electron. An example of such fiducial cut can be seen in 
Fig.[S]for a given electron momentum bin. The solid line in 
the top plot shows the boundary of the fiducial region for 
the central momentum in that bin. Only electron events 
inside the curve (blue area) were used in the analysis. This 
curve was determined by selecting the flat high-efficiency 
areas in the e -sliced 4> e distributions. The histograms on 
the bottom of Fig. show examples of such <j) e distribu- 
tions at two values of e = 23° ± 0.5° and 29° ± 0.5°. One 
sees a central, uniform area, flanked by two fringes. The 
highlighted area in the center indicates the selected fidu- 
cial range. In addition, a set of 9 e versus p e cuts was used 
to eliminate the areas with low detection efficiency due to 
problematic time-of- flight counters, photomultiplier tubes 
in Cherenkov counters, or drift chamber areas. 

4.2.2 Pion fiducial cuts 

The fiducial cuts for pions depend on the angles 6 V , tp^ 
and the momentum p^. The pion momentum is scanned in 
100 MeV steps from 0.3 to 1.7 GeV. The uniform detector 
efficiency region was determined by selecting a flat high- 
efficiency </>„■ region in each ^-sliced momentum bin, and 
the bad TOF counters and the inefficient DC areas were 
excluded by additional software cuts (the same procedure 
as was applied to electrons). Figure[7]shows an example for 




Fig. 5. (color online). Velocity /3 versus momentum for ir + 
candidates using experimental data (top) and the GSIM Monte 
Carlo simulation (bottom). The solid curves are ±1.5cr P cut 
lines used to select n + from positron (/3=1 band) and pro- 
ton (f3 <0.8 band) backgrounds. (The data are issued from 
a skimmed file which selected one electron and at least one 
positively charged particle and one neutral particle). 



the fiducial cuts for pions. The low-efficiency DC regions 
(between the black solid lines) and the bad TOF paddles 
(between red solid lines) are removed in both experimental 
(top) and simulated (bottom) data as part of the fiducial 
cuts. 

4.3 Kinematic corrections 

Due to effects that are not included in the reconstruc- 
tion software (deviations of the magnetic field from per- 
fect toroidal symmetry, misalignment of the tracking sys- 
tem,...), we have to apply some empirical corrections to 
the measured angles and momenta of both electrons and 
pions. For electrons, the kinematic corrections are applied 
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Fig. 6. (color online). Example of electron fiducial cuts for 
the electron momentum bin (p e = 1.437 ±.025 GeV) in Sector 
2. See the detailed explanation in the text. 




P K [GeV] 



Fig. 7. (color online). Pion polar angular distribution as a 
function of momentum in Sector 3. The low detector response 
areas are removed by empirical cuts for experimental (top) 
and simulated data (bottom). The black thin solid curves are 
fiducial cuts based on DC inefficiencies and the red thick solid 
curves are cuts for bad TOF counters. 



using the elastic ep — » e'p' process for which the kinemat- 
ics is over-constrained. The goal is to correct the three- 
momentum of the electron so as to minimize the con- 
straints due to the equations of conservation of energy and 
momentum. The same procedure is applied to the 7r + 's 
three-momentum using our reaction ep — > e'n + n under 
study, minimizing the deviation of the missing mass peak 
position from the neutron mass. The same correction fac- 
tors are used for all events having the same kinematics. In 
this way we keep the spatial resolution of the drift cham- 
ber systems and multiple scattering effects and the miss- 
ing mass resolution approaches its intrinsic limitations. 
The corrections were most sizable (w 5%) for the pion 
momentum. They resulted in an improved missing mass 
resolution, from 23 to 35 MeV depending on kinematics. 
The corrections were most sizable for the high-momentum 
and forward-angle pions at high W which are of interest in 
this experiment. We then applied additional ad-hoc smear- 



ing factors for the tracking and timing resolutions to the 
Monte Carlo so that they match the experimental data. 

5 Monte Carlo simulation 

In order to calculate the CLAS acceptance for ep — > e'ir + n, 
we simulated electron and pion tracks using the CLAS 
GEANT3-based Monte Carlo Package GSIM. For system- 
atic checks, we used two Monte Carlo event generators. 
Our approach is that by comparing the results of simu- 
lations carried out with two very different event gener- 
ators, a conservative and reliable estimation of system- 
atic effects, such as finite bin size effects, is obtained. The 
first event generator, GENEV (see Ref. [35] for the origi- 
nal publication dedicated to photoproduction processes), 
generates events for various exclusive meson electropro- 
duction reactions for proton and neutron targets (%, w, 
p , and (j>), including their decay, radiative effects, and 
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resonant and non-resonant multi-pion production, with re- 
alistic kinematic distributions. GENEV uses cross section 
tables based on existing photoproduction data and extrap- 
olates to electroproduction by introducing a virtual pho- 
ton flux factor (r) and the electromagnetic form factors. 
Radiative effects, based on the Mo and Tsai formula [25] . 
are part of this event generator as an option. Although the 
formula is exact only for elastic e-p scattering, it can be 
used as a first approximation to simulate the radiative tail 
and to estimate bin migration effects in our pion produc- 
tion process, as will be discussed in Sec. 15.21 The second 
event generator, FSGEN [30], distributes events according 
to the ep — > e'ir + n phase space. 

Electrons and positive pions were generated under the "el- 
6" experimental conditions. Events were processed through 
GSIM. As already mentioned, additional ad-hoc smear- 
ing factors for the tracking and timing resolutions are 
applied after GSIM so that they match the experimen- 
tal data. The low-efficiency regions in the drift chambers 
and problematic TOF channels were removed during this 
procedure. Acceptance and radiative corrections were cal- 
culated for the same kinematic bins as were used for the 
yield extraction as shown in Table [1] Figure [5] shows the 
binning in Q 2 and ig applied in this analysis. However, 
some bins will be dropped at some later stage in the anal- 
ysis, in particular due to very low acceptances (see follow- 
ing subsection). Our cross sections will be defined at the 
(xb,Q 2 , —t) values given by the geometrical center of the 
three-dimensional bins. To account for non-linear varia- 
tions of the cross section within a bin, a correction to our 
cross sections is determined by fitting with a simple ad- 
hoc three- variable function the simultaneous (xb, Q 2 , — in- 
dependence of our cross sections. This correction comes 
out at the level of a couple of percent in average. 



5.1 Acceptance corrections 

We related the experimental yields to the cross sections us- 
ing the acceptance, including the efficiency of the detector. 
The acceptance factor (Acc) compensates for various ef- 
fects, such as the geometric coverage of the detector, hard- 
ware and software inefficiencies, and resolution from track 
reconstruction. We generated approximately 850 million 
events, taking radiative effects into account. This results 
in a statistical uncertainty for the acceptance determina- 
tion of less than 5% for most bins, which is much lower 
than the systematic uncertainty that we have estimated 
(see Sec.©. 

We define the acceptance as a function of the kinematic 
variables, 



A < n 2 f A* "i N REC [x B ,Q\-t,<l>l) 



where N RE is the number of reconstructed particles and 
^rad N ls t ne number of generated particles in each kine- 
matic bin (the meaning of the subscript rad. will be- 
come clear in the next section) . The kinematic variables in 
jyGEN re f er t the generated values so that bin migration 
effects are taken into account in the definition of our accep- 
tance. The acceptances are in general between 1 and 9%. 
Figure [S] shows examples of acceptances, determined with 
the GENEV+GSIM packag function of the angle 

<j>* at a given Q 2 for various xb and t bins. Bins with an 
acceptance below 0.2% were dropped. For the integration 
over the 4>* angle, in order to obtain our three-fold cross 
sections, we fitted the acceptance-corrected 0* distribu- 
tions, so that any hole in the 0* distribution would be 
replaced by its fit value. 




0.2 0.4 0.6 



Fig. 8. (color online). Kinematic coverage and binning (red 
boxes) as a function of is and Q 2 (integrated over all other 
variables) for the experimental data. Only events with W > 
2 GeV are shown. 



5.2 Radiative correction 

Our goal is to extract the so-called Born cross section 
(tree-level) of the p(e, e'n + )n which can thus be compared 
to models. However, we measure a process which is accom- 
panied by higher order radiative effects. Our experimen- 
tal cross section must therefore be corrected. Radiative 
corrections are of two types: "virtual" corrections where 
there is no change in the final state of the p(e, e'n + )n pro- 
cess and "real" ones where there is in addition one (or 
several) Bremsstrahlung photon(s) in the final state. Such 
real Bremsstrahlung photons can originate either from the 
primary hard scattering at the level of the target proton 
(internal radiation) or from the interaction of the scattered 
or the initial electron with the various material layers of 
the CLAS detector that it crosses (external radiation). 

We have dealt with these corrections in two steps. The 
effects of the radiation of hard photons (for instance, the 
loss of events due to the application of a cut on the neu- 
tron missing mass) are taken into account by the Monte 
Carlo acceptance calculation described in the previous sec- 
tion. Indeed, as mentioned earlier, the GENEV code has 
the option to generate radiative photons according to the 
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Fig. 10. (color online). The simulated neutron missing mass 
distributions for two W bins with AW = 100 MeV at W = 
1.95 (left) and 2.15 (right) GeV integrated over <j)%, cos#J, and 
Q 2 . Normalized yields are shown with (solid red) and without 
(dashed blue) radiative effects. 



Fig. 9. (color online). Examples of acceptance as a function 
of </>* for various t and xb bins at Q 2 — 2.35 GeV 2 . The dips 
at <j>% = 0° and 180° are due to sector boundaries in CLAS. 



Mo and Tsai formula and the events in Eq. [5] were 

generated with this option turned on. Figure [TU] shows ex- 
amples of the simulated neutron missing mass with and 
without radiative effects in two W bins, obtained with 
the GENEV event generator and GSIM. Again, the Monte 
Carlo simulations were carried out with the same cut pro- 
cedures and conditions as used in the data analysis. 

Then, the correction due to soft photons and virtual 
corrections is determined by extracting the ratio between 
the number of events without radiative and with radiative 
effects at the level of GENEV for each three-dimensional 
kinematic bin. We therefore apply the following additional 
correction factor to our data: 



F rad { XBl Q , -t) N GEN {xB ^_ t y ( 3 ) 

As a check, these radiative-correction factors were also 
calculated with the EXCLURAD code [21] , which contains 
a complete description of all internal radiative effects in 
exclusive processes, but is currently valid only up to W = 
2 GeV. We compare the two different radiative-correction 
methods in a kinematic region where both methods are 
valid. Figure QT] shows the results for radiative-correction 
factors in the region W w 1.75 GeV and Q 2 m 3 GeV 2 as 
a function of cos 8* . 

The radiative correction factors from EXCLURAD are 
within ±20% of unity over the full cos 9* range (red solid 
points). The radiative corrections from GENEV+GSIM 
also fluctuate around 1.0 with a similar structure (blue 
open circles). The GENEV+GSIM error bars are due to 
Monte Carlo statistics (EXCLURAD is a theoretical code 



which has therefore no statistical uncertainty). The agree- 
ment between the two approaches is important because 
EXCLURAD is believed to be the most reliable of the two 
methods because it does not have the limitations of Mo 
and Tsai. Building on this reasonable agreement in this 
part of the phase space, we rely on the GENEV+GSIM 
radiative-correction factors for our data. In Sec. [3 we dis- 
cuss the systematic uncertainty associated with these ra- 
diative corrections. 




Fig. 11. (color online). Radiative-correction factors (RC) as 
a function of cos#* from EXCLURAD (red solid points) at 
W = 1.74 GeV, Q 2 = 3 GeV 2 , and cj>l = 112.5° and GENEV 
plus GSIM (blue open circle) at W « 1.75 GeV, Q 2 w 3 GeV 2 , 
and 80° < 4>l < 120°. 
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6 Background subtraction 

There are two main sources of background in our reaction. 
One consists of the mis-identification of pions with other 
positively charged particles (protons, kaons, positrons). 
This is particularly important for the pion-proton sepa- 
ration at high- momenta (p > 2 GeV), see Sec. 14.11 The 
other consists of multi-pion production. To subtract both 
backgrounds, we fit the neutron missing mass distribution 
bin by bin. We used many methods to fit these spectra: fit 
of only the background, fit of the signal plus background, 
with different functional forms both for the signal and the 
background, variation of the fitted range, etc... from which 
we extracted a systematic uncertainty (see Sec. [7]). 

Figure H2] (top) shows an example of a fit based on 
only the background, with an exponential plus a Gaussian. 
The former function was determined from simulations of 
the multi-pion spectra in the neutron missing mass region 
> 1.02 GeV. A comparison of the missing mass (MMx) 
spectrum is shown in the bottom plot of Fig. Q2] before 
(black squares) and after (red solid circles) background 
subtraction. In the range of the neutron missing mass 
cut, shown by the two vertical lines at 0.877 GeV and 
1.0245 GeV, the background is small, and the remaining 
radiative tail becomes visible after the background is sub- 
tracted. 



7 Systematic uncertainties 

Several sources of systematic uncertainty that can affect 
our measurements have been studied by changing various 
cuts and using different event generators. 

We varied the criteria used for the particle identifica- 
tion to provide more and less stringent particle selection 
simultaneously for experimental and GSIM data and then 
reran the complete analysis. The cuts on EC energy depo- 
sition and CC amplitude for the electron, as well as cuts 
on the TOF timing for the pion, have been varied. The EC 
sampling fraction cut was varied from ±3<tec to ±2ctec 
which led to a 5% uncertainty for electron identification. 
Changing the TOF (3 cut from ±2<jtof to ±2.5<7tof f° r 
pion identification gives a 1.7% uncertainty. The various 
cuts for channel identification such as fiducial, missing 
mass, and vertex cuts produced 3%, 1%, and 1.6% sys- 
tematic uncertainties, respectively. 

Acceptance and radiative corrections are the biggest 
sources of systematic errors. The systematic uncertainty 
from the acceptance is evaluated by comparing our results 
using the GENEV and FSGEN event generators. In the 
limit of infinitely large statistics and infinitely small bin 
size, our acceptances should be model-independent (up to 
the bin- migration effects). But these conditions are not 
reached here and we find differences between 2 and 8%. 
The systematic uncertainty for radiative corrections is es- 
timated similarly by comparing the radiative-correction 
factors (GENEV and EXCLURAD). We calculated the 
difference between the cross sections corrected for radia- 
tive effects using either GENEV-GEANT simulation or 
the W-expanded EXCLURAD (where EXCLURAD was 
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Fig. 12. (color online). Example of the missing mass peak plus 
background at Q 2 = 2.65 GeV 2 , -t = 1.15 GeV 2 , and x B = 
0.43. The top plot shows the fitted background distribution 
(hashed region). The bottom plot shows the neutron missing 
mass distribution before (black squares) and after (red solid 
points) background subtraction. 



linearly extrapolated to W > 2 GeV). An average 8% sys- 
tematic uncertainty was found. Acceptance and radiative 
corrections are actually correlated, but after a combined 
analysis we estimated an averaged range 4 — 12% total 
uncertainty for both of these effects together. 

Concerning the background subtraction procedure un- 
der the neutron missing mass (see Sec. |()]), we used vari- 
ous fitting functions (Gaussian plus exponential, Gaussian 
plus polynomial, exponential plus polynomial, etc.) and 
various fitting ranges. These various fitting functions and 
ranges eventually produced small differences and we esti- 
mated a 3% systematic uncertainty associated with this 
procedure. 

To take into account the model-dependency of our bin- 
centering correction (see Sec. l4.3p . we also introduce an er- 
ror equal to the correction factor itself which is, we recall, 
at the level of a couple of percent in average. 

These latter systematic uncertainties were determined 
for each bin. Concerning overall scale uncertainties, the 
target length and density have a 1% systematic uncer- 
tainty and the integrated charge uncertainty is estimated 
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at 2%. The background from the target cell was sub- 
tracted based on the empty-target runs and amounted to 
0.6±0.2% of our e'n + X events. The total systematic un- 
certainties, averaged over all bins, is then approximately 
12%. Table [5] summarizes the main systematic uncertain- 
ties in this analysis averaged over all the accessible kine- 
matic bins seen in Fig. [8] 

Table 2. Average systematic uncertainties for the differential 
cross sections. 



gion W > 2 GeV. We have extracted the differential cross 
sections as a function of several variables (t, Q 2 , and W or 
Xb)- The angle 0* is always integrated over in the follow- 
ing. The extraction of the interference cross sections gtt 
and (Jtl is the subject of an ongoing analysis and will be 
presented in a future article. The error bars on all cross 
sections include both statistical and systematic uncertain- 
ties added in quadrature. All values of our cross sections 
and uncertainties can be found on the CLAS database web 
page: |http://clasweb.jlab.org/cgi-bin/clasdb/db.cgi 



Source 


Criterion 


Estimated 
contribution 


Type 


point-to-point 




e~ PID 


i~ i ViLl 1 1 '1111 ^1 11 CLVj L L\J 11 

cut in EC 
(3asF — > 2o"sf) 


5% 


e — fiducial cut 


fiducial volume change 
(10% reduced) 


2.5% 


7T+ PID 


j3 resolution change 
(2<ttof - > 2.5otof) 


1.7% 


7r + fiducial cut 


width (10% reduced) 


3.5% 


Missing 
cut 


neutron missing 

rnflcq rPQollltlOn 

(3<tmmx — y 3.5<tmmx) 


1% 


Vertex cut 


2-vertex width 
(5% reduced) 


1.6% 


Acceptance 

Radiative 

corrections 


GENEV vs FSGEN 
GENEV vs EXCLURAD 


4-12% 


Background 
subtraction 


various fit functions 
exponential, gaussian 
and high order polynomials 


3% 


Bin-centering 
effect 


toy model 


2-4% 


Type 


overall scale/normalization 




LH2 target 


density /length 


1% 


Luminosity 


integrated charge 


2% 


Total 




9-14% 



8 Results and Discussion 

In this section, we present our results for the cross sec- 
tions of the p(e, eV + )n reaction in the invariant mass re- 



8.1 da /dt as a function of t 

Fig.rnHshows the differential cross section da/dt as a func- 
tion of t for different (xb, Q 2 ) bins. We define the reduced 
differential cross section: 



da 
~dl 



1 



d :i c 



r dQ 2 dxBdt 



(4) 



where the virtual photon flux factor [32] has been fac- 
tored out. We have included in Fig. [T3] the JLab Hall C 
data, which cover only the very small t domain. We note 
that there is generally reasonable agreement between the 
results of the two experiments. However, care must be 
taken in comparing the Hall C and Hall B measurements 
as the central (t, Q 2 , and W or xb, e) values do not ex- 
actly match each other. For instance, the most important 
discrepancy seems to appear in the bin (xb, Q 2 )=(0.49, 
3.35) where the Hall C measurement was carried out at 
e=0.45 [7| while ours corresponds to e=0.58 (the xb and 
Q 2 values being almost similar). According to the value 
of ctl relative to or, the Hall C cross section should then 
be renormalized: by a factor of 1.58/1.45~10% (if a^ax 
which the Hall C separated data [SUB] indicate, although at 
a slightly different kinematics) to a factor 0.58/0.45«30% 
(if ah dominates over ot which the Laget model predicts). 
For better visualization, which is also relevant for the com- 
parison with the models, we also show Fig. [14] which con- 
centrates on the low \t\ range of Fig. [TBI 

The da/dt cross sections fall as |i| increases, with some 
flattening at large \t\, which is a feature that is also ob- 
served in photoproduction [T1I20] . For several bins, for in- 
stance (x B , Q 2 )=(0.31, 1.75) or (0.37, 2.05), we notice a 
structure in da/dt for \t\ w 0.5 GeV 2 . The origin of this 
dip is not known. We note that the JLab Hall C exper- 
iment [7] also measured such a structure in da/dt (see 
Fig. 13 in Ref. [7 for bin (W, Q 2 )=(1.8, 2.16)). 

We first compare our data to calculations using hadronic 
degrees of freedom. The first one with which we will com- 
pare our data is the Laget model [33] based on Reggeized 



7r + and p + meson exchanges in the t— channel 



The 



hadronic coupling constants entering the calculation are 
all well-known or well-constrained, and the main free pa- 
rameters are the mass scales of the electromagnetic form 
factors at the photon-meson vertices. 

If one considers only standard, monopole, Q 2 -dependent 
form factors, one obtains much steeper t-slopes than the 
data. An agreement with the data can be recovered by 
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introducing a form factor mass scale that also depends on 
t according to the prescription of Ref . [33] . This form fac- 
tor accounts phenomenologically for the shrinking in size 
of the nucleon system as t increases. The size of the ef- 
fect is quantitatively the same as in the p(e, e'ui)p channel 
(see Fig. 1 of Ref. [53]), which is dominated by pion ex- 
change in the same energy domain as in our study. The 
results of this calculation with (Q 2 , independent meson 
electromagnetic form factors are shown, for daL/dt and 
da/dt = dar/dt + eddL/dt, in Figs. [T3l and [T4l by the red 
curves. The Laget model gives a qualitative description of 
the data, with respect to the overall normalization at low 
t and the xb~, Q 2 - and t- dependencies. We recall that this 
model already gives a good description of the photopro- 
duction data (SLAC, JLab) and of the HERMES electro- 
production data, and that the form factor mass scale [33] 
has not been adjusted to fit our data. 

In the framework of this model, daL/dt dominates at 
low \t\, while dax/dt takes over around |t| w 2 GeV 2 , 
this value slightly varying from one (Q 2 , xb) bin to an- 
other. This dominance of ol at low \t\ is a consequence 
of the i-channel 7r + -exchange (pion pole). At larger \t\, 
the p + meson exchange, which contributes mostly to the 
transverse part of the cross section, begins to dominate. 
The Laget Regge model, in addition to i-channel meson 
exchanges, also contains u-channel baryon exchanges. It 
thus exhibits an increase of the cross section in some (Q 2 , 
Xb) bins at the largest |i|-values, corresponding to low-|u| 
values. We have additional data at larger \t\ (lower 
that are currently under analysis. 

The second model with which we compare our data 
is the "hybrid" two-component hadron-parton model pro- 
posed in Refs. [551136) . Like in the Laget model, it is based 
on the exchange of the n + and p + Regge trajectories in the 
t-channel. However, the model complements these hadron- 
like interaction types, which dominate in photoproduc- 
tion and low Q 2 electroproduction, by a direct interaction 
of virtual photons with partons at high values of Q 2 fol- 
lowed by string (quark) fragmentation into 7r + n. The par- 
tonic part of the production mechanism is described by a 
"deep inelastic" -like electroproduction mechanism where 
the quark knockout reaction 7*0 — > q is followed by the 
fragmentation process of the Lund type. The transverse 
response is then treated as the exclusive limit of the semi- 
inclusive reaction p(e, e'ir + )X. Figures [13] and [14] show 
the results of this model compared to our data where very 
good agreement is found. This calculation was also found 
to give a good description of the L/T-separated Hall C 
and unseparated HERMES data [551155] . 

The third model that we wish to discuss, the GK model, 
is based purely on partonic degrees of freedom and is based 
on the handbag GPD formalism. In this model daL/dt is 
also mostly generated by the pion pole, similar to the two 
previous models. There are, however, a couple of differ- 
ences in the treatment of this pion pole in the GK calcula- 
tion. For instance, the Laget model has an intrinsic energy 
dependence. It is "Reggeized" , so the t-channel propaga- 
tor is proportional to s""^, where a n (t) is the pion Regge 
trajectory. In addition, it uses a (Q 2 , ^-dependent elec- 



tromagnetic form factor. These two features change the 
s-, xb~, and t- dependencies of the pion pole with respect 
to the GK treatment. Indeed, in the latter case, the t- 
channel pion propagator is proportional to l/(t — m 2 ), so 
it has no s-dependence, and the hadronic form factor at 
the ttNN vertex is only ^-dependent. 

Figure [15] shows the results of the GK calculation (in 
blue) for daL/dt and da/dt. We recall that the GK model 
is applicable only for small values of ~t/Q 2 . Outside this 
regime, higher-twist contributions that are not taken into 
account in the GK handbag formalism are expected. The 
GK model describes qualitatively our law-t unseparated 
cross sections over our whole (xb, Q 2 ) domain. This is 
remarkable since the GK model was optimized for higher- 
energy kinematics (HERMES) and no further adjustments 
were made for the present CLAS kinematics. We see that 
d<Jh / dt has a non- negligible contribution only in the very 
low \t\ domain and only for a few (xb, Q 2 ) bins, in par- 
ticular at the lowest xb and the largest Q 2 values (for 
instance, the (x B , Q 2 ) bins (0.25, 1.75) and (0.31, 2.35)). 
This is in line with the observation that at HERMES kine- 
matics, i.e. at lower xb and larger Q 2 values, the longitu- 
dinal part of the cross section dominates in the GK model 
at low \t\. For the larger xb values, one sees that the dom- 
inance of dohjdt at low \t\ is not at all systematic in the 
GK calculation. The ratio of daL/dt to da/dt strongly de- 
pends on xb ■ Specifically, it decreases as xb increases and 
at £#=0.49, daL/dt is only a few percent of da/dt, even 
at the lowest t values. This is a notable difference from 
the Laget Regge model for instance. 
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Fig. 13. (color online). Differential cross sections da/dt [/ib/GeV 2 ] integrated over <f>% for various (Q 2 , Xb) bins. The blue solid 
points are the present work. The error bars (outer error) on all cross sections include both statistical (inner error) and systematic 
uncertainties added in quadrature. The black open squares (da/dt) [7] and open stars (daL/dt) [5] are JLab Hall C data. The 
red thick solid (da/dt), and dashed (daL/dt) curves are the calculations from the Laget model |33| with (Q 2 , independent 
form factors at the photon-meson vertex. The black thin solid (da/dt) and dashed (daL/dt) curves are the calculations from 
the Kaskulov et al. model |35| . This model does not provide the calculation at (xb, Q 2 ) = (0.43, 2.35), (0.49, 2.95) due to the 
kinematic limit in the model. 
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Fig. 14. (color online). Same as Fig. ll3l except with an expanded low \t~t m in \ scale, where -tmin is the minimum kinematically 
possible four-momentum transfer. In addition, the black solid triangles [B] show the JLab Hall C extracted doL/dt data. 
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In particular, one can remark in Fig. 1151 where we dis- 
play in two (x B , Q 2 ) bins ((0.31, 1.75) and (0.37, 2.35)) 
the longitudinal part of the cross section as extracted from 
Hall C [5], that the longitudinal part of the GK calcu- 
lation is not in good agreement with the experimental 
data. This can be attributed to the way the pion pole 
and/or the pion-nucleon form factor, which are the main 
contributors to the longitudinal part of the cross section, 
are modeled in the GK approach. A Reggeization (like in 
the Laget model) or a change in the pion-nucleon form 
factor parametrization could possibly enhance the pion 
pole contribution at JLab kinematics and provide better 
agreement with our data (without damaging the agree- 
ment with the HERMES data) [37]. We recall that the 
GK model for which the GPD parameters were fitted to 
the low xb HERMES data, was simply extrapolated to 
the kinematics of the present article without any opti- 
mization and thus the present disagreement observed in 
dohjdt should not be considered as definitive. 

In the GK model, the transverse part of the cross sec- 
tion is due to transversity GPDs. In Fig. [T5j the GK cal- 
culation predicts that the transverse part of the cross sec- 
tion dominates essentially everywhere in our kinematic do- 
main. Although the GK L/T ratio probably needs to be 
adjusted as we just discussed, the GK calculation opens 
the original and exciting perspective to access transversity 
GPDs through exclusive 7r + electroproduction. 

Finally, at the kinematics of our experiment, in spite 
of our W > 2 GeV cut, it cannot be excluded that nu- 
cleoli resonances contribute. In Ref. [35], Kaskulov and 
Mosel identify these high-lying resonances with partonic 
excitations in the spirit of the resonance-parton duality 
hypothesis and invoke the continuity in going from an in- 
clusive final deep inelastic state to exclusive pion produc- 
tion. During this transition one expects that the inclusion 
of resonance excitations enhances the transverse response 
while leaving the longitudinal strength originating in the 
t-channel meson exchanges intact. Thus, in this work, the 
t-channel exchange part of the production amplitude is 
again described by the exchange of the Regge trajectories 
(ir + , p + and af) to which it is added a nucleon resonance 
component that is described via a dual connection be- 
tween the resonance and partonic deep inelastic processes. 
The parameters of this model have been tuned using the 
forward JLab Hall C data. Figure [15] shows the results of 
this calculation with our data and a reasonable agreement 
is found. 

The four models that we just discussed, although they 
give a reasonable description of the unseparated cross sec- 
tions, display rather different L/T ratios. The precise mea- 
surement of this ratio as a function of xb , Q 2 and t appears 
thus as essential to clarify the situation. For instance, in 
order to validate and/or tune the GK approach, it would 
be interesting to study the Q 2 -dependence at fixed xb and 
t of the longitudinal and transverse cross sections. They 
should approach, as Q 2 increases, a 1/Q 6 and 1/Q 8 scaling 
behavior respectively, as mentioned in the introduction of 
this article. In contrast, the Laget Regge model, for which 
xb is not a "natural" variable (it is rather W) should not 



predict such a Q 2 scaling at fixed xb- Although we are 
probably very far from such an asymptotic regime, the 
measurement of the Q 2 -dependence in the transition re- 
gion accessible with the upcoming JLab 12-GeV upgrade 
should provide some strong constraints and in particu- 
lar some checks on the way the higher-twist corrections 
are treated in the GK model. Such a program is already 
planned at JLab [55] . 

8.2 da/dt as a function of Q 2 at fixed t 

Figures [T75] and ITTlshow the differential cross section da/dt 
as a function of Q 2 at fixed xb for various t values. In 
Fig- HI °ur data arc fitted with a 1/Q" function and are 
compared to the GK model. We recall that, at asymptot- 
ically large Q 2 , the handbag mechanism predicts a dom- 
inance of ox which should scale as 1/Q 6 at fixed t and 
xb- The resulting exponents n of our fit indicates a flat- 
ter Q 2 dependence than 1/Q 6 . At the relatively low Q 2 
range accessed in this experiment, higher-twist effects are 
expected to contribute and hence the leading-twist 1/Q 6 
dependence of ctl is no longer expected. We note that such 
higher-twist contributions are part of the GK calculation 
and the GK model also does not show this scaling be- 
havior at the present Q 2 values. Although the GK model 
tends to underestimate the normalization of our data, its 
Q 2 dependence agrees reasonably well with our data. 

In Fig. [T71 we compare our data to the Laget [33] and 
the Kaskulov et al. 35,36 models. The Laget calcula- 
tion gives a reasonable description of the data although 
it seems to have a slightly steeper Q 2 -dependence than 
our data (particularly in the £s=0.37 bin). We note that 
in the xs=0.43 bin, our data seem to display a structure 
(dip) for Q 2 values between 3 and 4 GeV 2 , which is cer- 
tainly intriguing. We have at this stage no particular ex- 
planation for this. We just observe that the "hybrid" two- 
component hadron-parton model of Refs. 35,36 displays 
apparently also such structure which should therefore be 
further investigated. 

8.3 da/dt as a function of W at fixed 

Figure [15] shows our scaled cross sections, s 1 da/dt, as a 
function of W for four Q 2 values and four bins in cos#*: 
-0.01 ± 0.16, 0.27 ±0.1, 0.42 ± 0.05 and 0.53 ± 0.06. The 
lever arm in W is limited. At 9* = 90°, where the scal- 
ing behavior is expected to set in most quickly, we have 
only 2 or 3 data points in W, depending on the Q 2 bin. It 
is therefore difficult to draw precise conclusions at this 
stage for the VF-dependence at fixed Q 2 . Nevertheless, 
with these limited (but unique) data, one can say that, 
at 9* = 90°, except for the 3 data points at Q 2 =2.35 
GeV 2 , the VF-dependence of s 7 da/dt does not appear to 
be constant. We also display in Fig. [JS] the result of the 
Laget model. It gives, within a factor two, a general de- 
scription of these large-angle data. The W^- dependence of 
our data is very similar to the energy dependence that 
was observed in photoproduction [2]. In the same energy 
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range as covered by the present study, real-photon data 
exhibit strong deviations from scaling. Within the Laget 
model, these deviations are accounted for by the coupling 
between the nir + and the pN channels [3D]. The JLab 12- 
GeV upgrade will allow us to increase the coverage in W 
and check whether the hints of oscillations that we observe 
remain in the virtual-photon sector. 



9 Summary 

Wc have measured the cross sections of exclusive electro- 
production of 7r + mesons from protons as a function of 
-t = 0.1 - 5.3 GeV 2 , x B = 0.16 - 0.58, and Q 2 = 1.6 
-4.5 GeV 2 . We have compared our differential cross sec- 
tions to four recent calculations based on hadronic and 
partonic degrees of freedom. The four models give a qual- 
itative description of the overall strength and of the t-, Q 2 - 
and Xb~ dependencies of our unseparated cross sections. 
There is an obvious need for L-T separated cross sections 
in order to distinguish between the several approaches. 
These separations will be possible with the upcoming JLab 



12-GeV upgrade. In particular, if the handbag approach 
can accomodate the data, the p(e,e'ir + )n process offers 
the outstanding potential to access transversity GPDs. 
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